Comparison of a multiplex PCR with DNA barcoding for identification of container breeding mosquito species.

BACKGROUND: Identification of mosquitoes greatly relies on morphological specification. Since some species cannot be distinguished reliably by morphological methods, it is important to incorporate molecular techniques into the diagnostic pipeline. DNA barcoding using Sanger sequencing is currently widely used for identification of mosquito species. However, this method does not allow detection of multiple species in one sample, which would be important when analysing mosquito eggs. Detection of container breeding Aedes is typically performed by collecting eggs using ovitraps. These traps consist of a black container filled with water and a wooden spatula inserted for oviposition support. Aedes mosquitoes of different species might lay single or multiple eggs on the spatula. In contrast to Sanger sequencing of specific polymerase chain reaction (PCR) products, multiplex PCR protocols targeting specific species of interest can be of advantage for detection of multiple species in the same sample. METHODS: For this purpose, we adapted a previously published PCR protocol for simultaneous detection of four different Aedes species that are relevant for Austrian monitoring programmes, as they can be found in ovitraps: Aedes albopictus, Aedes japonicus, Aedes koreicus, and Aedes geniculatus. For evaluation of the multiplex PCR protocol, we analysed 2271 ovitrap mosquito samples from the years 2021 and 2022, which were collected within the scope of an Austrian nationwide monitoring programme. We compared the results of the multiplex PCR to the results of DNA barcoding. RESULTS: Of 2271 samples, the multiplex PCR could identify 1990 samples, while species determination using DNA barcoding of the mitochondrial cytochrome c oxidase subunit I gene was possible in 1722 samples. The multiplex PCR showed a mixture of different species in 47 samples, which could not be detected with DNA barcoding. CONCLUSIONS: In conclusion, identification of Aedes species in ovitrap samples was more successful when using the multiplex PCR protocol as opposed to the DNA barcoding protocol. Additionally, the multiplex PCR allowed us to detect multiple species in the same sample, while those species might have been missed when using DNA barcoding with Sanger sequencing alone. Therefore, we propose that the multiplex PCR protocol is highly suitable and of great advantage when analysing mosquito eggs from ovitraps.

A high heterozygosity genome assembly of Aedes albopictus enables the discovery of the association of PGANT3 with blood-feeding behavior.

The Asian tiger mosquito, Aedes albopictus, is a global invasive species, notorious for its role in transmitting dangerous human arboviruses such as dengue and Chikungunya. Although hematophagous behavior is repulsive, it is an effective strategy for mosquitoes like Aedes albopictus to transmit viruses, posing a significant risk to human health. However, the fragmented nature of the Ae. albopictus genome assembly has been a significant challenge, hindering in-depth biological and genetic studies of this mosquito. In this research, we have harnessed a variety of technologies and implemented a novel strategy to create a significantly improved genome assembly for Ae. albopictus, designated as AealbF3. This assembly boasts a completeness rate of up to 98.1%, and the duplication rate has been minimized to 1.2%. Furthermore, the fragmented contigs or scaffolds of AealbF3 have been organized into three distinct chromosomes, an arrangement corroborated through syntenic plot analysis, which compared the genetic structure of Ae. albopictus with that of Ae. aegypti. Additionally, the study has revealed a phylogenetic relationship suggesting that the PGANT3 gene is implicated in the hematophagous behavior of Ae. albopictus. This involvement was preliminarily substantiated through RNA interference (RNAi) techniques and behavioral experiment. In summary, the AealbF3 genome assembly will facilitate new biological insights and intervention strategies for combating this formidable vector of disease. The innovative assembly process employed in this study could also serve as a valuable template for the assembly of genomes in other insects characterized by high levels of heterozygosity.

The transcriptional response in mosquitoes distinguishes between fungi and bacteria but not Gram types.

Mosquitoes are prolific vectors of human pathogens, therefore a clear and accurate understanding of the organization of their antimicrobial defenses is crucial for informing the development of transmission control strategies. The canonical infection response in insects, as described in the insect model Drosophila melanogaster, is pathogen type-dependent, with distinct stereotypical responses to Gram-negative bacteria and Gram-positive bacteria/fungi mediated by the activation of the Imd and Toll pathways, respectively. To determine whether this pathogen-specific discrimination is shared by mosquitoes, we used RNAseq to capture the genome-wide transcriptional response of Aedes aegypti and Anopheles gambiae (s.l.) to systemic infection with Gram-negative bacteria, Gram-positive bacteria, yeasts, and filamentous fungi, as well as challenge with heat-killed Gram-negative, Gram-positive, and fungal pathogens. From the resulting data, we found that Ae. aegypti and An. gambiae both mount a core response to all categories of infection, and this response is highly conserved between the two species with respect to both function and orthology. When we compared the transcriptomes of mosquitoes infected with different types of bacteria, we observed that the intensity of the transcriptional response was correlated with both the virulence and growth rate of the infecting pathogen. Exhaustive comparisons of the transcriptomes of Gram-negative-challenged versus Gram-positive-challenged mosquitoes yielded no difference in either species. In Ae. aegypti, however, we identified transcriptional signatures specific to bacterial infection and to fungal infection. The bacterial infection response was dominated by the expression of defensins and cecropins, while the fungal infection response included the disproportionate upregulation of an uncharacterized family of glycine-rich proteins. These signatures were also observed in Ae. aegypti challenged with heat-killed bacteria and fungi, indicating that this species can discriminate between molecular patterns that are specific to bacteria and to fungi.

Mutations and intron polymorphisms in voltage-gated sodium channel genes of different geographic populations of Culex pipiens pallens/Culex pipiens quinquefasciatus in China.

BACKGROUND: Culex pipiens pallens and Culex pipiens quinquefasciatus are the dominant species of Culex mosquitoes in China and important disease vectors. Long-term use of insecticides can cause mutations in the voltage-gated sodium channel (vgsc) gene of mosquitoes, but little is known about the current status and evolutionary origins of vgsc gene in different geographic populations. Therefore, this study aimed to determine the current status of vgsc genes in Cx. p. pallens and Cx. p. quinquefasciatus in China and to investigate the evolutionary inheritance of neighboring downstream introns of the vgsc gene to determine the impact of insecticides on long-term evolution. METHODS: Sampling was conducted from July to September 2021 in representative habitats of 22 provincial-level administrative divisions in China. Genomic DNA was extracted from 1308 mosquitoes, the IIS6 fragment of the vgsc gene on the nerve cell membrane was amplified using polymerase chain reaction, and the sequence was used to evaluate allele frequency and knockdown resistance (kdr) frequency. MEGA 11 was used to construct neighbor-joining (NJ) tree. PopART was used to build a TCS network. RESULTS: There were 6 alleles and 6 genotypes at the L1014 locus, which included the wild-type alleles TTA/L and CTA/L and the mutant alleles TTT/F, TTC/F, TCT/S and TCA/S. The geographic populations with a kdr frequency less than 20.00% were mainly concentrated in the regions north of 38° N, and the geographic populations with a kdr frequency greater than 80.00% were concentrated in the regions south of 30° N. kdr frequency increased with decreasing latitude. And within the same latitude, the frequency of kdr in large cities is relatively high. Mutations were correlated with the number of introns. The mutant allele TCA/S has only one intron, the mutant allele TTT/F has three introns, and the wild-type allele TTA/L has 17 introns. CONCLUSIONS: Cx. p. pallens and Cx. p. quinquefasciatus have developed resistance to insecticides in most regions of China. The neighboring downstream introns of the vgsc gene gradually decreased to one intron with the mutation of the vgsc gene. Mutations may originate from multiple mutation events rather than from a single origin, and populations lacking mutations may be genetically isolated.

Molecular and entomological surveillance of malaria vectors in urban and rural communities of Benguela Province, Angola.

BACKGROUND: Malaria is a major public health problem in Angola, with Anopheles gambiae sensu lato (s.l.) and An. funestus s.l. being the primary vectors. This study aimed to clarify the information gaps concerning local Anopheles mosquito populations. Our objectives were to assess their abundance, geographical dispersion, and blood-feeding patterns. We also investigated their insecticide resistance. Molecular methods were used to identify sibling species, determine the origin of blood meals, measure Plasmodium falciparum infection rates, and detect the presence of knockdown resistance (kdr) mutations. METHODS: Adult mosquitoes were collected indoors using CDC light traps from nine randomly selected households at two sentinel sites with distinct ecological characteristics. The samples were collected from 1 February to 30 June 2022. Anopheles mosquitoes were morphologically identified and subjected to molecular identification. Unfed Anopheles females were tested for the presence of P. falciparum DNA in head and thorax, and engorged females were screened for the source of the blood meals. Additionally, members of An. gambiae complex were genotyped for the presence of the L1014F and L1014S kdr mutations. RESULTS: In total, 2226 adult mosquitoes were collected, including 733 Anopheles females. Molecular identification revealed the presence of Anopheles coluzzii, An. gambiae senso stricto (s.s.), An. arabiensis, and An. funestus s.s. Notably, there was the first record of An. coluzzii/An. gambiae s.s. hybrid and An. vaneedeni in Benguela Province. Plasmodium falciparum infection rates for An. coluzzii at the urban sentinel site and An. funestus s.s. at the rural site were 23.1% and 5.7%, respectively. The L1014F kdr mutation was discovered in both resistant and susceptible An. coluzzii mosquitoes, while the L1014S mutation was detected in An. gambiae s.s. for the first time in Benguela Province. No kdr mutations were found in An. arabiensis. CONCLUSIONS: This study provides valuable insights into the molecular characteristics of malaria vectors from the province of Benguela, emphasising the need for continuous surveillance of local Anopheles populations regarding the establishment of both kdr mutations for tailoring vector control interventions.

Extensive variation and strain-specificity in dengue virus susceptibility among African Aedes aegypti populations.

African populations of the mosquito Aedes aegypti are usually considered less susceptible to infection by human-pathogenic flaviviruses than globally invasive populations found outside Africa. Although this contrast has been well documented for Zika virus (ZIKV), it is unclear to what extent it is true for dengue virus (DENV), the most prevalent flavivirus of humans. Addressing this question is complicated by substantial genetic diversity among DENV strains, most notably in the form of four genetic types (DENV1 to DENV4), that can lead to genetically specific interactions with mosquito populations. Here, we carried out a survey of DENV susceptibility using a panel of seven field-derived Ae. aegypti colonies from across the African range of the species and a colony from Guadeloupe, French West Indies as non-African reference. We found considerable variation in the ability of African Ae. aegypti populations to acquire and replicate a panel of six DENV strains spanning the four DENV types. Although African Ae. aegypti populations were generally less susceptible than the reference non-African population from Guadeloupe, in several instances some African populations were equally or more susceptible than the Guadeloupe population. Moreover, the relative level of susceptibility between African mosquito populations depended on the DENV strain, indicating genetically specific interactions. We conclude that unlike ZIKV susceptibility, there is no clear-cut dichotomy in DENV susceptibility between African and non-African Ae. aegypti. DENV susceptibility of African Ae. aegypti populations is highly heterogeneous and largely governed by the specific pairing of mosquito population and DENV strain.

The updated genome of the Hungarian population of Aedes koreicus.

Vector-borne diseases pose a potential risk to human and animal welfare, and understanding their spread requires genomic resources. The mosquito Aedes koreicus is an emerging vector that has been introduced into Europe more than 15 years ago but only a low quality, fragmented genome was available. In this study, we carried out additional sequencing and assembled and characterized the genome of the species to provide a background for understanding its evolution and biology. The updated genome was 1.1 Gbp long and consisted of 6099 contigs with an N50 value of 329,610 bp and a BUSCO score of 84%. We identified 22,580 genes that could be functionally annotated and paid particular attention to the identification of potential insecticide resistance genes. The assessment of the orthology of the genes indicates a high turnover at the terminal branches of the species tree of mosquitoes with complete genomes, which could contribute to the adaptation and evolutionary success of the species. These results could form the basis for numerous downstream analyzes to develop targets for the control of mosquito populations.

A systematic review: is Anopheles vagus a species complex?

BACKGROUND: Anopheles vagus (subgenus Cellia) has been identified as a vector for malaria, filariasis, and Japanese encephalitis in Asia. Sporozoites of Plasmodium falciparum and Plasmodium vivax have been found in this zoophilic mosquito in Asia and Indonesia. This study systematically reviews publications regarding An. vagus species, variation, bio-ecology, and malaria transmission in various localities in Asia, especially Indonesia, to determine whether the current data support An. vagus as a species complex. METHODS: The databases Pubmed, Scopus, Europe PMC, and Proquest were searched to identify information regarding the morphology, karyotypes, polytene chromosome, cross-mating, ecology, and molecular identification of An. vagus was then evaluated to determine whether there were possible species complexes. RESULTS: Of the 1326 articles identified, 15 studies were considered for synthesis. The Anopheles spp. samples for this study came from Asia. Eleven studies used morphology to identify An. vagus, with singular studies using each of karyotype identification, chromosomal polytene identification, and cross-breeding experiments. Ten studies used molecular techniques to identify Anopheles spp., including An. vagus. Most studies discovered morphological variations of An. vagus either in the same or different areas and ecological settings. In this review, the members of An. vagus sensu lato grouped based on morphology (An. vagus, An. vagus vagus, An. vagus limosus, and An. limosus), karyotyping (form A and B), and molecular (An. vagus genotype A and B, An. vagus AN4 and AN5). Genetic analysis revealed a high conservation of the ITS2 fragment among members except for the An. vagus genotype B, which was, in fact, Anopheles sundaicus. This review also identified that An. vagus limosus and An. vagus vagus were nearly identical to the ITS2 sequence. CONCLUSION: Literature review studies revealed that An. vagus is conspecific despite the distinct morphological characteristic of An. vagus and An. limosus. Further information using another barcoding tool, such as mitochondrial COI and ND6 and experimental cross-mating between the An. vagus and An. limosus may provide additional evidence for the status of An. vagus as a species complex.

Structural and functional comparisons of salivary α-glucosidases from the mosquito vectors Aedes aegypti, Anopheles gambiae, and Culex quinquefasciatus.

Mosquito vectors of medical importance both blood and sugar feed, and their saliva contains bioactive molecules that aid in both processes. Although it has been shown that the salivary glands of several mosquito species exhibit α-glucosidase activities, the specific enzymes responsible for sugar digestion remain understudied. We therefore expressed and purified three recombinant salivary α-glucosidases from the mosquito vectors Aedes aegypti, Anopheles gambiae, and Culex quinquefasciatus and compared their functions and structures. We found that all three enzymes were expressed in the salivary glands of their respective vectors and were secreted into the saliva. The proteins, as well as mosquito salivary gland extracts, exhibited α-glucosidase activity, and the recombinant enzymes displayed preference for sucrose compared to p-nitrophenyl-α-D-glucopyranoside. Finally, we solved the crystal structure of the Ae. aegypti α-glucosidase bound to two calcium ions at a 2.3 Ångstrom resolution. Molecular docking suggested that the Ae. aegypti α-glucosidase preferred di- or polysaccharides compared to monosaccharides, consistent with enzymatic activity assays. Comparing structural models between the three species revealed a high degree of similarity, suggesting similar functional properties. We conclude that the α-glucosidases studied herein are important enzymes for sugar digestion in three mosquito species.

Successful introgression of wMel Wolbachia into Aedes aegypti populations in Fiji, Vanuatu and Kiribati.

Pacific Island countries have experienced periodic dengue, chikungunya and Zika outbreaks for decades. The prevention and control of these mosquito-borne diseases rely heavily on control of Aedes aegypti mosquitoes, which in most settings are the primary vector. Introgression of the intracellular bacterium Wolbachia pipientis (wMel strain) into Ae. aegypti populations reduces their vector competence and consequently lowers dengue incidence in the human population. Here we describe successful area-wide deployments of wMel-infected Ae. aegypti in Suva, Lautoka, Nadi (Fiji), Port Vila (Vanuatu) and South Tarawa (Kiribati). With community support, weekly releases of wMel-infected Ae. aegypti mosquitoes for between 2 to 5 months resulted in wMel introgression in nearly all locations. Long term monitoring confirmed a high, self-sustaining prevalence of wMel infecting mosquitoes in almost all deployment areas. Measurement of public health outcomes were disrupted by the Covid19 pandemic but are expected to emerge in the coming years.

Ecdysone-controlled nuclear receptor ERR regulates metabolic homeostasis in the disease vector mosquito Aedes aegypti.

Hematophagous mosquitoes require vertebrate blood for their reproductive cycles, making them effective vectors for transmitting dangerous human diseases. Thus, high-intensity metabolism is needed to support reproductive events of female mosquitoes. However, the regulatory mechanism linking metabolism and reproduction in mosquitoes remains largely unclear. In this study, we found that the expression of estrogen-related receptor (ERR), a nuclear receptor, is activated by the direct binding of 20-hydroxyecdysone (20E) and ecdysone receptor (EcR) to the ecdysone response element (EcRE) in the ERR promoter region during the gonadotropic cycle of Aedes aegypti (named AaERR). RNA interference (RNAi) of AaERR in female mosquitoes led to delayed development of ovaries. mRNA abundance of genes encoding key enzymes involved in carbohydrate metabolism (CM)-glucose-6-phosphate isomerase (GPI) and pyruvate kinase (PYK)-was significantly decreased in AaERR knockdown mosquitoes, while the levels of metabolites, such as glycogen, glucose, and trehalose, were elevated. The expression of fatty acid synthase (FAS) was notably downregulated, and lipid accumulation was reduced in response to AaERR depletion. Dual luciferase reporter assays and electrophoretic mobility shift assays (EMSA) determined that AaERR directly activated the expression of metabolic genes, such as GPI, PYK, and FAS, by binding to the corresponding AaERR-responsive motif in the promoter region of these genes. Our results have revealed an important role of AaERR in the regulation of metabolism during mosquito reproduction and offer a novel target for mosquito control.

Population genetic structure and wing geometric morphometrics of the filarial vector Armigeres subalbatus (Diptera: Culicidae) in Thailand.

Armigeres subalbatus (Diptera: Culicidae) is a mosquito species of significant medical and veterinary importance. It is widely distributed across Southeast and East Asia and is commonly found throughout Thailand. This study assessed the genetic diversity and population structure of Ar. subalbatus in Thailand using the cytochrome c oxidase subunit I (COI) gene sequences. Additionally, wing shape variations among these populations were examined using geometric morphometrics (GM). Our results demonstrated that the overall haplotype diversity (Hd) was 0.634, and the nucleotide diversity (π) was 0.0019. Significant negative values in neutrality tests (p < 0.05) indicate that the Ar. subalbatus populations in Thailand are undergoing a phase of expansion following a bottleneck event. The mismatch distribution test suggests that the populations may have started expanding approximately 16,678 years ago. Pairwise genetic differentiation among the 12 populations based on Fst revealed significant differences in 32 pairs (p < 0.05), with the degree of differentiation ranging from 0.000 to 0.419. The GM analysis of wing shape also indicated significant differences in nearly all pairs (p < 0.05), except for between populations from Nakhon Pathom and Samut Songkhram, and between those from Chiang Mai and Mae Hong Son, suggesting no significant difference due to their similar environmental settings. These findings enhance our understanding of the population structure and phenotypic adaptations of mosquito vectors, providing vital insights for the formulation of more efficacious vector control strategies.

The impact of agrochemical pollutant mixtures on the selection of insecticide resistance in the malaria vector Anopheles gambiae: insights from experimental evolution and transcriptomics.

BACKGROUND: There are several indications that pesticides used in agriculture contribute to the emergence and spread of resistance of mosquitoes to vector control insecticides. However, the impact of such an indirect selection pressure has rarely been quantified and the molecular mechanisms involved are still poorly characterized. In this context, experimental selection with different agrochemical mixtures was conducted in Anopheles gambiae. The multi-generational impact of agrochemicals on insecticide resistance was evaluated by phenotypic and molecular approaches. METHODS: Mosquito larvae were selected for 30 generations with three different agrochemical mixtures containing (i) insecticides, (ii) non-insecticides compounds, and (iii) both insecticide and non-insecticide compounds. Every five generations, the resistance of adults to deltamethrin and bendiocarb was monitored using bioassays. The frequencies of the kdr (L995F) and ace1 (G119S) target-site mutations were monitored every 10 generations. RNAseq was performed on all lines at generation 30 in order to identify gene transcription level variations and polymorphisms associated with each selection regime. RESULTS: Larval selection with agrochemical mixtures did not affect bendiocarb resistance and did not select for ace1 mutation. Contrastingly, an increased deltamethrin resistance was observed in the three selected lines. Such increased resistance was not majorly associated with the presence of kdr L995F mutation in selected lines. RNA-seq identified 63 candidate resistance genes over-transcribed in at least one selected line. These include genes coding for detoxification enzymes or cuticular proteins previously associated with insecticide resistance, and other genes potentially associated with chemical stress response. Combining an allele frequency filtering with a Bayesian FST-based genome scan allowed to identify genes under selection across multiple genomic loci, supporting a multigenic adaptive response to agrochemical mixtures. CONCLUSION: This study supports the role of agrochemical contaminants as a significant larval selection pressure favouring insecticide resistance in malaria vectors. Such selection pressures likely impact kdr mutations and detoxification enzymes, but also more generalist mechanisms such as cuticle resistance, which could potentially lead to cross-tolerance to unrelated insecticide compounds. Such indirect effect of global landscape pollution on mosquito resistance to public health insecticides deserves further attention since it can affect the nature and dynamics of resistance alleles circulating in malaria vectors and impact the efficacy of control vector strategies.

A genotyping array for the globally invasive vector mosquito, Aedes albopictus.

BACKGROUND: Although whole-genome sequencing (WGS) is the preferred genotyping method for most genomic analyses, limitations are often experienced when studying genomes characterized by a high percentage of repetitive elements, high linkage, and recombination deserts. The Asian tiger mosquito (Aedes albopictus), for example, has a genome comprising up to 72% repetitive elements, and therefore we set out to develop a single-nucleotide polymorphism (SNP) chip to be more cost-effective. Aedes albopictus is an invasive species originating from Southeast Asia that has recently spread around the world and is a vector for many human diseases. Developing an accessible genotyping platform is essential in advancing biological control methods and understanding the population dynamics of this pest species, with significant implications for public health. METHODS: We designed a SNP chip for Ae. albopictus (Aealbo chip) based on approximately 2.7 million SNPs identified using WGS data from 819 worldwide samples. We validated the chip using laboratory single-pair crosses, comparing technical replicates, and comparing genotypes of samples genotyped by WGS and the SNP chip. We then used the chip for a population genomic analysis of 237 samples from 28 sites in the native range to evaluate its usefulness in describing patterns of genomic variation and tracing the origins of invasions. RESULTS: Probes on the Aealbo chip targeted 175,396 SNPs in coding and non-coding regions across all three chromosomes, with a density of 102 SNPs per 1 Mb window, and at least one SNP in each of the 17,461 protein-coding genes. Overall, 70% of the probes captured the genetic variation. Segregation analysis found that 98% of the SNPs followed expectations of single-copy Mendelian genes. Comparisons with WGS indicated that sites with genotype disagreements were mostly heterozygotes at loci with WGS read depth < 20, while there was near complete agreement with WGS read depths > 20, indicating that the chip more accurately detects heterozygotes than low-coverage WGS. Sample sizes did not affect the accuracy of the SNP chip genotype calls. Ancestry analyses identified four to five genetic clusters in the native range with various levels of admixture. CONCLUSIONS: The Aealbo chip is highly accurate, is concordant with genotypes from WGS with high sequence coverage, and may be more accurate than low-coverage WGS.

The frequency of kdr and ace-1 alleles in Anopheles gambiae s.l. before and during indoor residual spraying (IRS) implementation and four years after IRS withdrawal in three districts in Atacora, Benin.

BACKGROUND: Indoor residual spraying (IRS) was first implemented in the Atacora department, Benin from 2011 to 2012 using bendiocarb (carbamate) followed by annual spraying with pirimiphos-methyl (organophosphate) from 2013 to 2018. Before and after IRS implementation in Atacora, standard pyrethroid insecticide-treated bed nets were the main method of vector control in the area. This study investigated the knockdown resistance (kdr) gene (L1014F) and the acetylcholinesterase (ace-1) gene (G119S), before and during IRS implementation, and 4-years after IRS withdrawal from Atacora. This was done to assess how changes in insecticide pressure from indoor residual spraying may have altered the genotypic resistance profile of Anopheles gambiae s.l. METHOD: Identification of sibling species of An. gambiae s.l. and detection of the L1014F mutation in the kdr gene and G119S mutation in ace-1 genes was done using molecular analysis. Allelic and genotypic frequencies were calculated and compared with each other before and during IRS implementation and 4 years after IRS withdrawal. The Hardy-Weinberg equilibrium and genetic differentiation within and between populations were assessed. RESULTS: Prevalence of the L1014F mutation in all geographic An. gambiae s.l. (An. gambiae s.s., Anopheles. coluzzii, Anopheles. arabiensis, and hybrids of "An. gambiae s.s. and An. coluzzii") populations increased from 69% before IRS to 87% and 90% during and after IRS. The G119S allele frequency during IRS (20%) was significantly higher than before IRS implementation (2%). Four years after IRS withdrawal, allele frequencies returned to similar levels as before IRS (3%). Four years after IRS withdrawal, the populations showed excess heterozygosity at the ace-1 gene and deficit heterozygosity at the kdr gene, whereas both genes had excess heterozygosity before and during IRS (FIS < 0). No genetic differentiation was observed within the populations. CONCLUSIONS: This study shows that the withdrawal of IRS with bendiocarb and pirimiphos-methyl may have slowed down the selection of individual mosquitoes with ace-1 resistance alleles in contrast to populations of An. gambiae s.l. with the L1014F resistance allele of the kdr gene. This may suggest that withdrawing the use of carbamates or organophosphates from IRS or rotating alternative insecticides with different modes of action may slow the development of ace-1 insecticide-resistance mutations. The increase in the prevalence of the L1014F mutation of the kdr gene in the population, despite the cessation of IRS, could be explained by the growing use of pyrethroids and DDT in agriculture and for other domestic use. More observational studies in countries where carbamates or organophosphates are still being used as public health insecticides may provide additional insights into these associations.

The population genetics of partial diapause, with applications to the aestivating malaria mosquito Anopheles coluzzii.

Diapause, a form of dormancy to delay or halt the reproductive development during unfavourable seasons, has evolved in many insect species. One example is aestivation, an adult-stage diapause enhancing malaria vectors' survival during the dry season (DS) and their re-establishment in the next rainy season (RS). This work develops a novel genetic approach to estimate the number or proportion of individuals undergoing diapause, as well as the breeding sizes of the two seasons, using signals from temporal allele frequency dynamics. Our modelling shows the magnitude of drift is dampened at early RS when previously aestivating individuals reappear. Aestivation severely biases the temporal effective population size ( N e $$ {N}_e $$ ), leading to overestimation of the DS breeding size by 1 / 1 - α 2 $$ 1/{\left(1-\alpha \right)}^2 $$ across 1 year, where α $$ \alpha $$ is the aestivating proportion. We find sampling breeding individuals in three consecutive seasons starting from an RS is sufficient for parameter estimation, and perform extensive simulations to verify our derivations. This method does not require sampling individuals in the dormant state, the biggest challenge in most studies. We illustrate the method by applying it to a published data set for Anopheles coluzzii mosquitoes from Thierola, Mali. Our method and the expected evolutionary implications are applicable to any species in which a fraction of the population diapauses for more than one generation, and are difficult or impossible to sample during that stage.

Blood-feeding patterns of Culex pipiens biotype pipiens and pipiens/molestus hybrids in relation to avian community composition in urban habitats.

BACKGROUND: Culex pipiens sensu stricto (s.s.) is considered the primary vector of Usutu virus and West Nile virus, and consists of two morphologically identical but behaviourally distinct biotypes (Cx. pipiens biotype pipiens and Cx. pipiens biotype molestus) and their hybrids. Both biotypes are expected to differ in their feeding behaviour, and pipiens/molestus hybrids are presumed to display intermediate feeding behaviour. However, the evidence for distinct feeding patterns is scarce, and to date no studies have related differences in feeding patterns to differences in host abundance. METHODS: Mosquitoes were collected using CO2-baited traps. We collected blood-engorged Cx. pipiens/torrentium specimens from 12 contrasting urban sites, namely six city parks and six residential areas. Blood engorged Cx. pipiens/torrentium mosquitoes were identified to the species and biotype/hybrid level via real-time polymerase chain reaction (PCR). We performed blood meal analysis via PCR and Sanger sequencing. Additionally, avian host communities were surveyed via vocal sounds and/or visual observation. RESULTS: We selected 64 blood-engorged Cx. pipiens/torrentium mosquitoes of which we successfully determined the host origin of 55 specimens. Of these, 38 belonged to biotype pipiens, 14 were pipiens/molestus hybrids and the identity of three specimens could not be determined. No blood-engorged biotype molestus or Cx. torrentium specimens were collected. We observed no differences in feeding patterns between biotype pipiens and pipiens/molestus hybrids across different habitats. Avian community composition differed between city parks and residential areas, whereas overall avian abundance did not differ between the two habitat types. CONCLUSIONS: Our results show the following: (1) Cx. pipiens s.s. feeding patterns did not differ between city parks and residential areas, regardless of whether individuals were identified as biotype pipiens or pipiens/molestus hybrids. (2) We detected differences in host availability between city parks and residential areas. (3) We show that in both urban habitat types, biotype pipiens and pipiens/molestus hybrids fed on both mammalian and avian hosts. This underscores the potential role in arbovirus transmission of biotype pipiens and pipiens/molestus hybrids.

Pyrethroid resistance status and co-occurrence of V1016G, F1534C and S989P mutations in the Aedes aegypti population from two dengue outbreak counties along the China-Myanmar border.

BACKGROUND: Over the past two decades, dengue fever (DF) has emerged as a significant arboviral disease in Yunnan province, China, particularly in the China-Myanmar border area. Aedes aegypti, an invasive mosquito species, plays a crucial role in transmitting the dengue virus to the local population. Insecticide-based vector control has been the primary tool employed to combat DF, but the current susceptibility status of Ae. aegypti to commonly used insecticides is unknown. Assessment of Ae. aegypti resistance to pyrethroid insecticides and an understanding of the underlying mechanisms of this resistance in the China-Myanmar border region is of significant strategic importance for effectively controlling the DF epidemic in the area. METHODS: Aedes aegypti larvae collected from Ruili and Gengma counties in Yunnan Province were reared to adults in the laboratory and tested for susceptibility to three pyrethroid insecticides (3.20% permethrin, 0.08% lambda-cyhalothrin and 0.20% deltamethrin) by the standard WHO susceptibility bioassay. Genotyping of mutations in the knockdown gene (kdr), namely S989P, V1016G and F1534C, that are responsible for resistance to pyrethroid insecticides was performed using allele-specific PCR methods. A possible association between the observed resistant phenotype and mutations in the voltage-gated sodium channel gene (VGSC) was also studied. RESULTS: Aedes aegypti mosquitoes collected from the two counties and reared in the laboratory were resistant to all of the pyrethroids tested, with the exception of Ae. aegypti from Gengma County, which showed sensitivity to 0.20% deltamethrin. The mortality rate of Ae. aegypti from Ruili county exposed to 3.20% permethrin did not differ significantly from that of Ae. aegypti from Gengma County (χ2 = 0.311, P = 0.577). By contrast, the mortality rate of Ae. aegypti from Ruili County exposed to 0.08% lambda-cyhalothrin and 0.20% deltamethrin, respectively, was significantly different from that of Ae. aegypti from Gengma. There was no significant difference in the observed KDT50 of Ae. aegypti from the two counties to various insecticides. Four mutation types and 12 genotypes were detected at three kdr mutation sites. Based on results from all tested Ae. aegypti, the V1016G mutation was the most prevalent kdr mutation (100% prevalence), followed by the S989P mutation (81.6%) and the F1534C mutation (78.9%). The constituent ratio of VGSC gene mutation types was significantly different in Ae. aegypti mosquitoes from Ruili and those Gengma. The triple mutant S989P + V1016G + F1534C was observed in 274 Ae. aegypti mosquitoes (60.8%), with the most common genotype being SP + GG + FC (31.4%). The prevalence of the F1534C mutation was significantly higher in resistant Ae. aegypti from Ruili (odds ratio [OR] 7.43; 95% confidence interval [CI] 1.71-32.29; P = 0.01) and Gengma (OR 9.29; 95% CI 3.38-25.50; P = 0.00) counties than in susceptible Ae. aegypti when exposed to 3.20% permethrin and 0.08% lambda-cyhalothrin, respectively. No significant association was observed in the triple mutation genotypes with the Ae. aegypti population exposed to 3.20% permethrin and 0.20% deltamethrin resistance (P > 0.05), except for Ae. aegypti from Gengma County when exposed to 0.08% lambda-cyhalothrin (OR 2.86; 95% CI 1.20-6.81; P = 0.02). CONCLUSIONS: Aedes aegypti from Ruili and Gengma counties have developed resistance to various pyrethroid insecticides. The occurrence of multiple mutant sites in VGSC strongly correlated with the high levels of resistance to pyrethroids in the Ae. aegypti populations, highlighting the need for alternative strategies to manage the spread of resistance. A region-specific control strategy for dengue vectors needs to be implemented in the future based on the status of insecticide resistance and kdr mutations.

Molecular monitoring of insecticide resistance in major disease vectors in Armenia.

BACKGROUND: Armenia is considered particularly vulnerable to life-threatening vector-borne diseases (VBDs) including malaria, West Nile virus disease and leishmaniasis. However, information relevant for the control of the vectors of these diseases, such as their insecticide resistance profile, is scarce. The present study was conducted to provide the first evidence on insecticide resistance mechanisms circulating in major mosquito and sand fly populations in Armenia. METHODS: Sampling sites were targeted based mainly on previous historical records of VBD occurrences in humans and vertebrate hosts. Initially, molecular species identification on the collected vector samples was performed. Subsequently, molecular diagnostic assays [polymerase chain reaction (PCR), Sanger sequencing, PCR-restriction fragment length polymorphism (RFLP), quantitative PCR (qPCR)] were performed to profile for major insecticide resistance mechanisms, i.e. target site insensitivity in voltage-gated sodium channel (vgsc) associated with pyrethroid resistance, acetylcholinesterase (ace-1) target site mutations linked to organophosphate (OP) and carbamate (CRB) resistance, chitin synthase (chs-1) target site mutations associated with diflubenzuron (DFB) resistance and gene amplification of carboxylesterases (CCEs) associated with resistance to the OP temephos. RESULTS: Anopheles mosquitoes were principally represented by Anopheles sacharovi, a well-known malaria vector in Armenia, which showed no signs of resistance mechanisms. Contrarily, the knockdown resistance (kdr) mutations V1016G and L1014F/C in the vgsc gene were detected in the arboviral mosquito vectors Aedes albopictus and Culex pipiens, respectively. The kdr mutation L1014S was also detected in the sand fly, vectors of leishmaniasis, Phlebotomus papatasi and P. tobbi, whereas no mutations were found in the remaining collected sand fly species, P. sergenti, P. perfiliewi and P. caucasicus. CONCLUSIONS: This is the first study to report on molecular mechanisms of insecticide resistance circulating in major mosquito and sand fly disease vectors in Armenia and highlights the need for the establishment of systematic resistance monitoring practices for the implementation of evidence-based control applications.